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^ \ Abstract 

^ ■ We propose a first principle computation of the equilibrium thermodynamics of simple 

\ fragile glasses starting from the two body interatomic potential. A replica formulation 

translates this problem into that of a gas of interacting molecules, each molecule being 
built of m atoms, and having a gyration radius (related to the cage size) which vanishes 
at zero temperature. 

We use a small cage expansion, valid at low temperatures, which allows to compute 
the cage size, the specific heat (which follows the Dulong and Petit law), and the config- 
Q . urational entropy. 
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1 Introduction 



Take a three dimensional classical system consisting of many particles, interacting through a 
short range potential with a repulsive core. Very often this system will undergo, upon cooling 
or upon compression, a solidification into an amorphous solid state - the glass state. The 
conditions required for observing this glass phase is the avoidance of crystallisation, which can 
always be obtained through a fast enough quench (the meaning of 'fast' depends very much 
of the type of system) Q. There also exist potentials which naturally present some kind of 
frustration with respect to the crystalline structure. Whether their actual stable state is a 
crystal or a glass is not known, but they are known to solidify into glass states, even when 
cooled slowly - such is the case for instance of binary mixtures of hard spheres, soft spheres, 
or Lennard- Jones particles with appropriately different radii. These have been studied a lot in 
recent numerical simulations [@, f| |J || |f. 

Our aim is to compute the equilibrium thermodynamic properties of this glass phase, using 
the statistical mechanical approach, namely starting from the microscopic Hamiltonian (an 
attempt to build up a non equilibrium thermodynamic phenomenology can be found in 0). We 
shall therefore assume that crystallisation has been avoided, and consider only the amorphous 
solid state. The general framework of our approach finds its roots in old ideas of Kauzman 
||, Adam and Gibbs ||, which received a boost when Kirkpatrick, Thirumalai and Wolynes 
underlined the analogy between structural glasses and some generalized spin glasses ITTfl . This 



framework should provide a good description of fragile glass-formers. These are the systems 
in which the increase of relaxation time upon decreasing the temperature is much faster than 
Arrhenius- often parametrized as a Vogel-Fulcher law, displaying a divergence of the relaxation 
time at a finite temperature 0. In this approach the glass transition, measured from dynamical 
effects, is supposed to be associated with an underlying thermodynamic transition at the 
Kauzman or Vogel-Fulcher temperature T K . This ideal glass transition is the one which should 
be observed on infinitely long time scales in fragile glass-formers 0. This transition is of an 
unusual type, since it presents two apparently contradictory features: 

1. The transition is continuous (second order) from the thermodynamical point of view: 
the internal energy is continuous, and the transition is signalled by a discontinuity of the 
specific heat which jumps from its liquid value above Tk to a value very close to that of 
a crystal phase below. 



2. The order parameter is discontinuous at the transition. 



In order to make this last statement precise we shall have to define an order parameter 
for the glass phase in the framework of equilibrium statistical mechanics, which involves some 
subtleties and will be addressed below. At this introductory stage let us take loosely as an 
order parameter the correlation in the positions of the particles at very large times. In the 
liquid there is no correlation. In the glass the positions are correlated in time. Clearly the 
order parameter jumps discontinuously between the liquid phase and the glass phase. The 



two properties above are indeed observed in generalized spin glasses ||12[| . The problem of the 
existence or not of a diverging correlation length is still an open one fT3|j . 

This analogy is suggestive, but it also hides some very basic differences, like the fact that 
spin glasses have quenched disorder while structural glasses do not. The recent discovery of 
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some generalized spin glass systems without quenched disorder |TJ], [15], [16| has given credit 
to the idea that this analogy is not fortuitous. The problem was to turn this general idea 
into a consistent computational scheme allowing for some quantitative predictions. Important 
steps in this direction were invented in |17|, |TJ| , which showed how useful it is to study several 



coupled copies of the same system in order to characterize properly the glass phase. In a 
previous preliminary study, we used some of these ideas to estimate the glass temperature, 

However the approximations we did were not adequate for 
Here we concentrate instead on the properties 



arriving from the liquid phase [T9[ 
the description of the low temperature phase, 
of the glass phase itself and we introduce approximations which are much more appropriate 
to describe its properties particularly at low temperatures. We are now able to construct 
analytical tools for doing computations in the glass phase and to test the results in numerical 
(and eventually real) experiments. A brief description of a part of the present work has 
appeared in [pO] . 

In the next section we shall present in more details the general physical picture underlying 
our approach. In sect. 3 we shall explain why and how one can characterize and study the 
glass phase using a replicated liquid. Sect. 4 derives the Hamiltonian of the molecular liquid, 
which is studied in the next two sections, first of all by a small cage expansion in sect. 5, 
then by a molecular HNC closure in sect. 6. In sect. 7 we present the results of these various 
approximations concerning the glass transition temperature and the thermodynamic quantities. 
Sect. 8 gives a list of some directions into which this work could be extended. Two appendices 
contain the derivation of the molecular HNC closure on one hand, and its expansion to second 
order in the small cage parameter on the other hand. 



2 The basic scenario 

In this section we want to present some of the general ideas which provide a background 
to our approach. These have to do with the existence of a configurational entropy, and the 
identification of the glass transition as a point where the configurational entropy vanishes. 
These ideas are presented in general, without special reference to a specific system. They can 
be derived in great details in some mean field spin glass models. Although the microscopic 
description of these models is somewhat remote from the actual glass problem which interests 
us, we have included for completeness a short summary of some of the results found in these 
systems. This will help to formulate the basic hypotheses of our approach. 

2.1 Configurational entropy 

We consider a system of iV particles moving in a volume V of a d-dimensional space, and 
interacting by some short range potential. These could be for instance hard spheres or Lennard- 
Jones particles. 

Let us introduce the free energy functional F[p] which depends on the local particle density 
p(x) and on the temperature. We suppose that at sufficiently low temperature this functional 
has an exponentially large number of minima More precisely, the number of free energy 
minima with free energy density / = F/N is supposed to be exponentially large in some region 
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Figure 1: Qualitative shape of the configurational entropy versus free energy. The whole curve 
depends on the temperature. The saddle point which dominates the partition function, for m 
constrained replicas, is the point /* such that the slope of the curve equals m/T (for the usual 
unreplicated system, m = 1). If the temperature is small enough the saddle point sticks to the 
minimum / = / m j n and the system is in its glass phase. 

of free energies, f min {T) < f < f ma x{T): 

Af(f,T,N)nexp(NS c (f,T)) . (1) 

Exactly at zero temperature these minima coincide with the mimima of the potential energy 
as function of the coordinates of the particles. The function S c is called the complexity or 
the configurational entropy (it is the contribution to the entropy coming from the existence of 
an exponentially large number of locally stable configurations). The number of local minima 
is supposed to vanish outside of the region f m i n (T) < f < f ma x(T), and the configurational 
entropy S c (f,T) is supposed to go to zero continuously at f m i n (T), as found in all existing 
models so far (see figjl]). 

Let us first discuss the properties of the system at thermal equilibrium: we thus consider the 
case where each configuration of the system is assigned a probability given by its Boltzmann 
weight. We label the free energy minima by an index a. To each of them we can associate 
a free energy F a and a free energy density f a = F a /N. In the low temperature region we 
suppose that the total free energy of the system ($) can be well approximated by the sum of 
the contributions to the free energy of each particular minimum: 

Z = exp(-/3iV$) ~ £exp(-/3iV/ a ). (2) 

a 

For large values of iV we can write 

/fmax 
dfexp[-N(J3f-S c (f,T))]. (3) 
^ in in 

We can thus use the saddle point method and approximate the integral with the integrand 
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evaluated at its maximum. We find that 



$ = mm$(f) = f*-TS e (f*,T), (4) 

where 

^(f) = f-TS c (f,T). (5) 

This formula is quite similar to the usual formula for the free energy ,i.e. / = minE(E—TS(E)), 
where S(E) is the entropy density as a function of the energy density (E). The main difference 
is the fact that the total entropy of the system has been decomposed into the contribution due 
to small fluctuations around a given configuration (this piece has been included into /), and 
the contribution due to the existence of a large number of locally stable configurations, the 
configurational entropy. 

Calling /* the value of / which minimize $(/), we have two possibilities: 

• The minimum lies inside the interval and it can be found as the solution of the equation 
j3 = dS c /df. In this case we have 

$ = f*-TS:, S* c =S c (f*,T). (6) 

The system may stay in one of the many possible minima. The number of accessible 
minima is exp(NS*) . The entropy of the system is thus the sum of the entropy of a 
typical minimum and of S*, which is the contribution to the entropy coming from the 
exponentially large number of metastable configurations. 

• The minimum is at the extreme value of the range of variability of /: it sticks at /* = f m i n 
and the total free energy is $ = f m in- in this case the contribution of the configurational 
entropy to the free energy is zero. The different states which contribute to the free energy 
have a difference in free energy density which is of order iV _1 (a difference in total free 
energy of order 1). This situation is often encountered in spin glasses, both in usual 



cases of spin glasses with quenched disorder |22], ^3J, and also in some spin glass systems 
without quenched disorder |T3|, jn| Pf] , 



One aim of the theory of glasses at equilibrium could be to demonstrate from first principles 
the existence of a configurational entropy function such as depicted in fig. |], and to compute 
it. This is difficult to achieve. For instance Kepler's conjecture, a simple zero temperature 
statement saying that there is no denser packing of hard spheres in three dimensions than 



the fee lattice, has resisted a proof for more than three centuries ||24|1 . Here we shall take a 
more modest starting point: we shall assume the existence of the local minima and of the 
configurational entropy function with the general properties depicted above, and within this 
assumption we shall show how to compute (approximately but with a rather good accuracy, 
and one which can be improved systematically) the various properties of the system, including 
the configurational entropy function itself. 

2.2 Mean field situation 

So far, the only systems for which the above program could be carried out in all details are 
some type of mean field spin glasses with a discontinuous jump of the order parameter at the 



transition 12, 11, 25, BE 27, 28 . 
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Although we will not need all the ingredients that have been found in these other problems, 
it is useful to recall some of them; later on we will mention how this picture might be modified 
in a realistic - non mean field - system. The configurational entropy function is convex, and 
previous work indicates that it depends smoothly on the temperature, the main effect of a 
temperature change being a global shift of the free energies. Starting from high temperatures, 
we thus encounter the following temperature regions (we use here the language of liquids and 
glasses) . 

• For T > To the free energy functional is dominated by the uniform density solution, 
p(x) = p (there may exist close to To other minima [^, |26|], but their total contribution 
has a higher free energy than the uniform solution). The system is obviously in the fluid 
phase. 

• In the region where To > T > T%, the minimum of the function $(/) is within the 
interval [f min (T), f m ax(T)}. Therefore the system can stay in one of many different states. 
The entropy of the equilibrium system receives a contribution from the configurational 
entropy, which is non zero. A very surprising result, found in all generalized mean field 
spin glasses with discontinuous transition so far, is that the total free energy of the system 
including the configurational entropy contribution, $(/*), is equal to the free energy of the 



fluid solution with uniform p |T7], [L8[ . This result has not received a general explanation 
beyond the simple idea of the transition at To being a fragmentation of accessible phase 
space into many separated pockets, the total volume of which should be continuous at 
To- Although the thermodynamics is still given by the usual expressions of the liquid 
phase and the free energy is analytic at T D , below this temperature the system, at each 
given moment of time, may stay in one of the exponentially large number of minima. 

• In the region where T < Tk the saddle point of $ sticks at its minimum and the free 
energy is dominated by the contribution of a few minima having the lowest possible value 
fmin{T). Here the free energy is no more the analytic continuation of the free energy in 
the fluid phase. A phase transition is present at Tk and the specific heat is discontinuous 
here. 

The intermediate phase To > T > Tk is particularly interesting. In the mean field systems, 
an exact solution of the Langevin dynamics indicates a dynamical phase transition at To, the 
system being trapped in some states with a free energy which is extensively higher than that 



of the equilibrium state |3(|. For the realistic finite dimensional problems which we want to 
study, the situation is much less clear, but one can speculate that the system will equilibrate 
in this regime, very slowly fll |. The time to jump from one minimum to another minimum is 



quite large: it is an activated process which is controlled by the height of the barriers which 
separate the different minima. The correlation time will become very large below To and for 
this reason To is called the dynamical transition point. The correlation time (which should be 
proportional to the viscosity) diverges only at the true thermodynamic transition temperature, 
sometimes called the ideal glass temperature Tk (see fig.0). The precise form of this divergence 
is not well understood. It is natural to suppose that one should get a divergence of the form 
exp(A/(T — Tk) u ) for an appropriate value of v, but a reliable analytic computation of v is 
lacking JTTl |3T| . Experiments can often be fitted by this law with various values of v, including 
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Figure 2: Relaxation time vs temperature in discontinuous spin-glasses. The right hand curve 
is the mean-field prediction, which gives a dynamical transition at a temperature To above 
the static transition temperature Tk- The left curve is a conjecture on the behaviour in finite 
dimensional systems: activated processes smear the dynamic transition. The relaxation time 
diverges only at the static temperature Tk, but becomes experimentally large already around 
the glass temperature T g . 

the Vogel-Fulcher fit with v — 1. The equilibrium configurational entropy is different from 
zero (and it is a number of order 1) when the temperature is smaller than To, it decreases 
with the temperature and it vanishes linearly at T = Tk- At this temperature the entropy of a 
single minimum becomes equal to the total entropy and the contribution of the configurational 
entropy to the total free energy vanishes. Therefore the total entropy and the internal energy 
are continuous at the transition. 

2.3 Relationship to experiments 

The above scenario is appealing in that it puts into a unified framework a number of experi- 
mental facts on glasses, as well as some general theoretical ideas. 

Experimentally, the system falls out of equilibrium when its relaxation time becomes larger 
than the experimental time. The 'glass transition temperature', defined conventionally as the 
temperature where the typical relaxation time reaches a value of order one hour, falls some- 
where between Tk and To- By considering slower and slower quenches, one can equilibrate the 
system at lower temperatures. However in this scenario there exists an underlying thermody- 
namic transition at the temperature TV, which is the ideal glass transition temperature. This 
temperature is also the one where the viscosity would diverge in the Vogel-Fulcher type fitting 
of the viscosity versus temperature. Clearly it also corresponds to the Kauzman temperature: 
the excess entropy of the supercooled liquid with respect to the crystal is basically equal to 
the configurational entropy, which vanishes precisely at Tk- The experimental fact that the 
Kauzman temperature and the Vogel-Fulcher one are close to each other has been noted many 
times, and is also found in the Adam-Gibbs relation |J. 

The dynamical temperature To also receives a natural interpretation. In mean field, there- 
fore neglecting activated processes, the relaxation time diverges with a power law at To, and 
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the autocorrelation function develops an infinitely long plateau. This slowing down is described 
precisely by the mode coupling theory ||10|| . In the mean field approximation the height of the 
barriers separating the different minima is infinite and the temperature Tjj is sharply defined 
as the point where the correlation time diverges. In the real world activated process (which are 
neglected in the mean field approximation and consequently in the mode coupling theory) have 
the effect of producing a finite (but large) correlation time also at and below To (the precise 
meaning of the dynamical temperature beyond mean field approximation is not so clear -see 
; probably the best definition is that To is the temperature where the mode-coupling theory 
predicts a transition). Therefore one expects that the mode coupling description will give good 
results in the region largely above Tjj, a fact that has been checked accurately in experiments 
P^l and numerical simulations [|33| . 

A last point which is predicted within the basic scenario, and has been checked numerically, 
is a specific type of aging and modification of the fluctuation-dissipation relation. The aging 
behaviour, which has been seen many years ago already in some polymeric glasses [53], can be 



studied in details in spin glasses [35]. These studies, initiated by the work of Cugliandolo and 



Kurchan ]30|], have led to some well defined generalisation of the basic equilibrium properties 



such as time translation invariance and fluctuation-dissipation theorem (FDT). This generali- 
sation is not limited to the narrow scope of some special mean field spin-glasses, but seems to 
provide a general description of glassy dynamics in many systems, including structural glasses. 
The modification of the fluctuation-dissipation relation can be measured directly, although the 
experiments are not simple. On the other hand, numerical simulations for a binary mixture of 
soft spheres [@ or Lennard- Jones particles |4| have found exactly the non-trivial modification 
which is predicted by the general scenario, providing therefore a confirmation of its validity at 
least on their (limited) time scales. 



3 A static order parameter for the glass phase 

In this section we wish to explain the general strategy for describing and computing properties 
of an amorphous solid state. We are particularly interested in systems with many metastable 
states, having a non zero configurational entropy. We shall explain the general strategy trying 
to keep away as much as possible from any specific model, the more precise formulation for our 
problem will be given in the next section. Let us consider a system of N particles, interacting 
by a two body potential with a Hamiltonian 

H= v{xi-Xj) (7) 

l<i<j<N 

where the particles move in a volume V of a d-dimensional space, and v is an arbitrary short 
range potential with a short range repulsion, like a 1/r 12 potential or a Lennard- Jones one. We 
shall take the thermodynamic limit N, V — > oo at fixed density p = N/V. For simplicity, we do 
not consider here the description of mixtures of different types of particles. The generalization 
to mixtures is necessary if one wants to compare more precisely to simulations, which are 
performed on mixtures in order to avoid crystallisation. This generalization, together with a 
detailed comparison, will be presented in a forthcoming paper fl40 ]. Some general background 
is provided by the review paper . 
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3.1 Time persistent correlations 

Before going to a purely static description of the order parameter, let us first discuss a dynami- 
cal one. At an atomic level one often tends to associate the glass transition with the divergence 
of the time scale on which a labelled particle can get out of its trap. While this is an intuitive 
picture, it is not possible to translate it into a good definition of the solid phase: because of the 
excitation and movements of vacancies and other defects, this individual trapping time scale is 
always finite, although it will increase exponentially when the temperature gets small. In order 



to get a proper definition of the solid, it has been proposed |4l], (42] to use a generalisation of 
the Edwards Anderson order parameter of the type: 

Q EA {p) = lim lira i £< e *G*<o)-. fc <*))) (8) 

jk 

where p is an arbitrary non zero wave vector, the order of magnitude of which is one over 
the typical interparticle distance. When the system is in the liquid phase the above order 
parameter is zero and when it is in the glass phase this order parameter is non zero (even in 
the presence of single particle diffusion). 

This definition would hold for the equilibrium dynamics, i.e. assuming that the system is in 
equilibrium at time t = 0. As we know the glass never reaches equilibrium and therefore it ages: 
correlations are not stationary in time. The proper generalization of the previous correlation 
taking into account the aging effect takes the slightly more complicated form (where the order 
of limits is crucial): 

Qea(p) — h m lim \im —Y(e ip - {x > itw) - Xk(tw+T)) ) (9) 

jk 

This gives a sensible dynamical definition of the glass phase. 
3.2 Correlations between two copies 

We would like a purely static description of the solid phase in the framework of equilibrium 
statistical mechanics, in a case where there are no Bragg peaks. As soon as we have a solid phase 
the translational symmetry is broken and the system can be in many states. For crystalline 
order these many states just differ from each other by rotations or translations which can be 
easily taken care of by appropriate boundary terms. In the glass case, in order to choose a 
state, one should first know the average position of each atom in the solid, which requires 
an infinite amount of information. Had we known this information, we could have added to 
the Hamiltonian an infinitesimal but extensive pinning field which attracts each particle to its 
equilibrium position, sending N to infinity first, before taking the limit of zero pinning field. 
This is the usual way of identifying the phase transition. 

In order to get around the problem of the description of the amorphous solid phase, a simple 
method has been developed in the spin glass context. Pictorially, one could say that although 
we do not know the conjugate field, the system itself knows it. The idea, borrowed from 
spin glass theory ||3, 44], is then to consider two copies of the system, with an infinitesimal 



extensive attraction. One then identifies the transition temperature from the fact that the two 



9 



replicas remain close to each other in the limit of vanishing coupling (having sent N to infinity 
first). 

In the case of glasses we can thus consider two identical systems of particles, {xj} and {%•}, 
with a total energy function: 

e = J2 ( v ( x i- x j) + v (yi-yj)) + e J2 w ( x i~yj) ( 10 ) 

where we have introduced a small attractive potential w(r) between the two systems. The 
precise shape of w is irrelevant, insofar as we shall be interested in the limit e — > 0, but its 
range should be of order or smaller than the typical interparticle distance. The order parameter 
is then the correlation function between the two systems: 




In the liquid phase this correlation function is identically equal to one, while it has a nontrivial 
structure in the glass phase, reminiscent of the pair correlation of a dense liquid, but with an 
extra peak around r ~ 0. Let us notice that we expect a discontinuous jump of this order 
parameter at the transition, in spite of its being second order in the thermodynamic sense. 
The existence of a non trivial order parameter is associated with the spontaneous breaking of 
a symmetry: For e = 0, with periodic boundary conditions, the system is symmetric under 
a global translation of the x particles with respect to the y particles. This symmetry is 
spontaneously broken in the low temperature phase, where the particles of each subsystem 
tend to sit in front of each other. One could equally use the Fourier transform of this cross- 
correlation, which then gives back, but in an equilibrium framework, the Edwards Anderson 
order parameter defined in (|). 



3.3 Thermodynamics below T#: replicas 



The previous method is a reasonable definition of an equilibrium order parameter which can be 
used in simulations or in analytic studies in order to identify the phase transition arriving from 
the liquid phase. However this technique can be improved in order to study quantitatively the 
glass phase itself. 

Let us assume that in the glass phase there exists a non zero configurational entropy, 
as introduced above. Clearly the knowledge of this configurational entropy as a function of 
free-energy and temperature, S c (f, T), will allow us to reconstruct all the interesting thermody- 
namic properties of the system. It has been realised by Monasson|T7| that the configurational 
entropy can be reconstructed from a study of an arbitrary number, m, of copies of the system, 
when they are constrained to be in the same state. As we will need to analytically continue 
the results in m, we shall call the copies 'replicas'. An alternative and related method is to 
introduce a real coupling of the system to another system which is thermalized [|18 1 ; this has 
been used recently in order to study the glass phase 0, |3 . The formulation which we present 



here is slightly different from, but equivalent to, that of pTfl . 

The basic idea is extremely simple. Instead of two copies of the system, let us consider m 
copies which are constrained to stay in the same minimum. We shall discuss below how one 
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can achieve this constraint, but let us first discuss the physics of this constrained system. Its 
partition function is: 

< jM df e -Nlmf3f-S c (f,T)] (12) 



J ra 



frt 



The dependence on the number m of replicas of the total free energy, 

$(m, T) = --^ log Z m « min(m / - TS c (f, Tj) , (13) 

allows to compute the configurational entropy S c (f, T) as a function of the free energy, using: 

d$(m,T) 



dm 
m 2 d<f)(m, T) 

T dm 

where 4>(m, T) is the free energy per particle: 



./ 

S c , (14) 



. $(m,T) , N 

m 

If the glass transition is due to the entropy crisis described in the previous section (and 
this is our main hypothesis), then the crucial quantity is the value of the slope Sq(T) of the 
configurational entropy at the lowest free energy: 

»oCT) = jf(fo(T)) (16) 

The usual glass transition is determined by TkSq{Tk) = 1. For the replicated and constrained 
system, the phase transition temperature depends on the number m of replicas and is 
determined by (see fig. [I]): 

T (m) So ( T (m)) = m ( 17 ) 

It is very natural to assume that sq(T) is a smooth function of temperature, going to a constant 
at zero temperature (we shall check this hypothesis self-consistently later). Then we see that, 
when m is continued analytically to real values, smaller than unity, one can have T( m ) < T K . 
The replicated and constrained system can thus be in the liquid phase for temperatures smaller 
than the glass transition temperature Tk = T^: it is then made up of molecules, each of which 
contains one atom of each replica, but these molecules are in a liquid state. The basic reason for 
this crucial fact is that for m < 1 the effective interaction potential (assuming for simplicity 
molecules of very small radius) is decreased from v(r) to mv(r), thus displacing the glass 
transition to lower temperatures. 

We are interested in the free energy in the glass phase, therefore in the region m = 1 and 
T < Tk- This free energy cannot be computed from that of the liquid with m = 1, T > Tk 
because of the phase transition at Tk- However we shall now show that one can deduce it 
from the free energy of the molecular fluid at m < 1. This molecular fluid with m < 1 has a 
transition to a glass state at the temperature T = < Tk- Inside the glass phase, thus for 
T < T^ m \ the free energy of the replicated and constrained system is given by the condition 

S c (f,T) = (18) 
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Figure 3: A sketch of the typical phase diagram in the temperature-m plane, for a system with 
m weakly coupled replicas. In the whole high temperature region above the full line the system 
is in a liquid phase. There are two liquid phases, above the horizontal line T = T K the various 
replicas are not correlated in the limit of the coupling e going to zero. On the contrary the 
liquid state at m < 1, in the region between the full line and the dashed line, is a molecular 
liquid where the various replicas form molecular bound states. The low temperature region 
below the full line (characterized by m = m*(T)) is the glass phase. In this glass phase, for 
a given temperature, the free energy per replica is m independent. Therefore one can deduce 
the free energy of the glass ( with m < 1 and T < T K ) from the knowledge of the free energy 
in the molecular liquid. 
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and it is independent on m. 

Let us now look at the phase diagram at a fixed temperature T < Tx, varying m (see fig. 
|^). The free energy per particle 4>(m, T) of the molecular liquid is an increasing function of m 
at small m, which reaches a maximum at a point m* < 1 where the glass transition takes place 
(obviously m* is the solution of: T Sq(T) = m*). As the free energy in the glass phase is m 
independent, the liquid free energy at the transition <p(m*,T) (which is equal to the glass free 
energy at the transition) is equal, for T < Tk, to the free energy <p(m = 1,T) of the glass at 
the temperature T. We have thus shown that the knowledge of the free energy of the molecular 
liquid, <p(m,T), allows to compute the free energy of the glass. 

These basic observations are at the heart of our strategy for computing properties of the 
glass phase. We shall write down more explicit formulas in our case below. We would like first 
to make three comments on this approach: 

• For T < Tk and m > m*, the free energy 4>(m) is constant and larger than the analytic 
continuation of the free energy 4>(m) of the molecular liquid. If one would have followed 
this molecular liquid in the region m > m*, one would have found that dip /dm < 0, 
predicting a negative configurational entropy. Instead, the glass transition occurs and 
the configurational entropy sticks to zero in the whole glass phase. The fact that the 
free energy in the glass phase is larger than the analytic continuation from the high 
temperature phase explains why the specific heat has a discontinuity downward when 
we decrease the temperature. This is in variance with what happens generally in other 
transitions (at least in the mean field approximation) where the free energy in the low 
temperature phase is smaller than the analytic continuation from the high temperature 
phase and the specific heat has a discontinuity upward when we decrease the temperature. 

• In practice in order to try to constrain the systems to be in the same state, one introduces 
some small attractive coupling, of order e, between the replicas. It is thus important to 
understand when this coupling leads to a molecular liquid. The phase diagram shown in 
fig|| can be conjectured from the following elementary study of the free energy, confirmed 
by exact computations of mean field discontinuous spin glasses [26, 18, [5|, 28]. There are 



a priori four possible cases. If the m replicas are in the same state, the free energy is 
$ = minf (m / — TS c (f, T)) — mim — l)e. If they are in different states, the free energy 
is $ = mmhiy-( / — TS c (f,T)). On top of this, the free energy minimum can either 
stick to fo (glass phase) or be at a value / larger than f (liquid). One just needs to 
find out which situation actually minimizes the free energy, for given values of m and T. 
The solution is displayed in fig|| showing that there is an intermediate molecular liquid 
phase at m < 1. 

The 'replicas' which we introduce here play a slightly different role compared to the ones 
used in disordered systems: there is no quenched disorder here, and no need to average 
a logarithm of the partition function. 'Replicas' are introduced to handle the problem 
of the absence of description of the amorphous state. We do not know of any other 
procedure to characterize an amorphous solid state in the framework of equilibrium 
statistical mechanics. There is no 'zero replica' limit, but there is, as in disordered 
systems, an analytic continuation in the number of replicas. We shall see that this 
continuation looks rather innocuous. 
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4 The replica approach to structural glasses: general 
formalism 



In this section we write down the formulas corresponding to the replica approach introduced in 
the previous section. We keep here to the case of simple glass formers consisting of N particles 
interacting by a pair potential v(r) in a space of dimension d. 

4.1 The partition function 

The usual partition function, used e.g. in the liquid phase, is 

1 r N 

z ^mJ Ifl^O e ~" H ( 19 ) 

We wish to study the transition to the glass phase through the onset of an off-diagonal cor- 
relation in replica space. We use m replicas and introduce the Hamiltonian of the replicated 
system: 

m 

H m= E Y,«-^)+ E in-!: o:j (20) 

l<i<j<JVa=l j • {I V} 

where W is an attractive interaction. The precise form of W is unimportant: it should be a 
short range attraction respecting the replica permutation symmetry, and its strength which 
will be sent to zero in the end. For instance one could take 

W{r\...,r m ) = e £ w(r a - r b ) (21) 

l<a<b<m 

with w(r) a smooth short range two body attraction. 
The partition function of the replicated system is 

i . N m 

nn^)e* (22) 

iv - J i=\a=\ 

The order parameter is the generalised cross correlation: 

p(r\...,0 = 1 £ < 8{xl -r*)...8{xZ -r™) > (23) 

jl—jm 

where the average is the Boltzmann-Gibbs average with the measure proportional to exp(—[3H m ). 

4.2 Molecular bound states 

At low enough temperature, we expect that the particles in the different replicas will stay close 
to each other due to the joint effect of the small inter-replica attraction and the intra-replica 
interactions: when the system is in the glass phase, the role of the attractive term W will be to 
insure that all replicas fall into the same glass state, so that the particles in different replicas 
stay at the same place, apart from some thermal fluctuations: A vanishingly small interaction 
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between replicas will give rise to a strong correlation. As the thermal fluctuations are relatively 
small throughout the solid phase (one can see this for instance from the Lindeman criterion), 
one can identify the molecules and relabel all the particles in the various replicas in such a 
way that the particle j in replica a always stays close to particle j in replica b. All the other 
relabelings are equivalent to this one, producing a global factor N\ m ~ l in the partition function. 

We therefore need to study a system of molecules, each of them consisting of m atoms (one 
atom from each replica). It is natural to write the partition function in terms of the variables 
Ti which describe the centers of masses of the molecules, and the relative coordinates uf, with 
xf = Ti + u\ and J2a u 1 = ; 

1 „ N N m N / \ 

z m = ^/nwnnwn ^e«?) 

iV • J i=l i=l o=l i=l V a ) 

exp [-(3 v{n - r, + < - u]) - P^Wtf, uf) J (24) 



5 The small cage expansion 

In order to transform these ideas into a tool for doing explicit computations of the thermody- 
namic properties of a glass we have to use an explicit method for computing the free energy as 
function of the temperature and m. As is usually the case, in the liquid phase exact analytic 
computations are not possible and we have to do some approximations. In this section we 
shall use the fact that the thermal fluctuations of the particles in the glass are small at low 
enough temperature: the size of the 'cage' seen by each particle is therefore small, allowing for 
a systematic expansion. What we will be describing here are the thermal fluctuations around 
the minimum of the potential of each particle, in the spirit of the Einstein model for vibrations 
of a crystal. 



5.1 Legendre transform 

We start from the replicated partition function Z m described in molecular coordinates in 
Assuming that the relative coordinates u" are small, we can expand W to leading order and 
write: 



1 „ N N m N I 

JnNniiwii ^(eo 

V • J i=l i=l o=l i=l \ a 



L — L L — L LL— ± l — ± \ 

exp L/3 £ v( ri - r, + < - «J) - ±- £ - n^) 2 ) (25) 

\ i<j,a a i a,b / 

In the end we are interested in the limit (1/a) — * 0. We would like first to define the size A of 
the molecular bound state, which is also a measure of the size of the cage seen by each atom 
in the glass, by: 

d(l/a) 2 * < fa 
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(d is the dimension, N is the number of particles). We Legendre transform the free energy 
4>(m,a) = —(T/m) \ogZ m , introducing the thermodynamic potential per particle ip{m, A): 

Mm, A) — Mm, a) + Td^—-^- — (27) 

2 a 

What we want to see is whether there exists a minimum of ip at a finite value of A. 

At low temperatures, this minimum should be at small A, and so we shall seek an expansion 
of ip in powers of A. It turns out that this can be found by an expansion of in powers of 
a, used as an intermediate bookkeeping in order to generate the low temperature expansion. 
This may look confusing since we are eventually going to send a to oo. However this method 
is nothing but a usual low temperature expansion in the presence of an infinitesimal breaking 
field. For instance if one wants to compute the low temperature expansion of the magnetization 
in a rf-dimensional Ising model in an infinitesimal positive magnetic field h, the main point is 
that the magnetisation is close to one. One can organise the expansion by studying first the 
case of a large magnetic field, performing the expansion in powers of exp(— 2h), and in the 
end letting h — > 0. A little thought shows that the intermediate -large h- expansion is just a 
bookkeeping device to keep the leading terms in the low temperature expansion. What we do 
here is exactly similar, the role of h being played by 1/a. 



5.2 Zeroth order term 

We use the equivalent form: 

1 r N m d d X- ( m \ 

z m = - } / n n ( rfd o n ex p -0 e - **) - v- - Xi)2 \ ■ (28) 



For a — * 0, the identity 



exp - X t f) c i^) d ' 2 n^ - XA (29) 



gives: 



Z ^ a)= V^) O ^y/n^exp -^E^-^) • (30) 

In this expression we recognise the integral over the X^s as the partition function Z Hq (T*) of 
the liquid at the effective temperature T* , defined by 

T* = T/m . (31) 

Therefore the free energy, at this leading order, can be written as: 

a 10/ \ d(l-m) 2ira d , . 1 . . . . 

m, a) = -±- log — log (m) log Z Uq (T ) 32 

2m m 2m mN 
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5.3 First order term 

In order to expand to next order, we start from the representation (|25|) and expand the inter- 
action term to quadratic order in the relative coordinates: 



7 



J nA/<n| m d 5£<)j exp (-pmY^viri-rj)] 

exp f - f E E« - - ^T^Mi - »v) - j- EK - «-) 2 

(The indices fi and z/, running from 1 to d, denote space directions). Notice that in order 
to carry this step, we need to assume that the interaction potential v(r) is smooth enough, 
excluding hard cores. To expand at small a we need the properties of the set of m random 
variables u a living on one site with measure P(u) cx 5(J2a ua ) exp(— (l/4a) J2ab( ua ~ u>> ) 2 )- It 
turns out that these are gaussian random variables with a first moment which vanishes and a 
second moment which is equal to: 



(«;«£>„= (S ab --)-5, u . (33) 



Expanding (j33f ) to first order in a we have: 

P 
2 



\ogZ m = log Zl - P - E E(K " u *n< - n]r)o(d,d u v(ri ~ r s ))* (34) 



4 <j aju> 

where the average (.)o is that for the u variables with the gaussian measure (^), and the average 
(.}* is over the center of mass positions r^, which are those of a liquid phase thermalized at the 
temperature T* = T/m. 

The free energy to first order is equal to: 

„ / n d(m—l) 1 „^ d(l — m) , 2n d , 1 , 

% a) = -\ log --a/3C+ -\ log — - — log m - — - log Z liq (T*) 35 

2m a 2m m 2m mJS 

where the constant C is proportional to the expectation value of the Laplacian of the potential, 
in the liquid phase at the temperature T*: 

C=~^Z(M«-*i)) m (36) 

£i III ■ f i -\ 

Differentiation with respect to - gives the size of the cage: 

_ 36 (1 — m). 2n ~ (1 — m) . . tnf _ s 

^^RT-\= ~n L da + a 2 (3C = - K ' dA 37 

o(l/ a) 2m 2 

Expanding this equation in perturbation theory in A we have: 

2(3m 3 C a2 , , 

a = mA- -j- -A 38 

d(m — 1) 
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The Legendre transform is then easily expanded to first order in A: 

Pip(m, A) = (3(f)(m, a) + d- 



2 a 



d(l — m) , . . ^ d(l — m) d , 1 , ^ ,_ Y , 

= log(2vrA) - /3mAC + >- - — log m - — - log Z Jiff T* 39 

2m 2m 2m miv 

This very simple expression gives the free energy as a function of the number of replicas, m, 
and the cage size A. We need to study it at m < 1, where we should maximise it with respect 
to A and m. The fact that we seek a maximum when m < 1 instead of the usual procedure of 
minimising the free energy is a well established fact of the replica method, appearing as soon 



as the number of replicas is less than 1 [22 



As a function of A , the thermodynamic potential ip has a maximum at: 

d(l — m) 1 d 1 

max = 2/3m 2 c = Pfd d rg*(r)Av(r) ^ ' 

where g* is the pair correlation of the liquid at the temperature T*. A study of the potential 
tfj(m, A max ), which equals 4>(m), as a function of m then allows to find all the thermodynamic 
properties which we seek, using the formulas of the previous section. This step and the results 
will be explained below in sect. [7], where we shall also compare the results to those of other 
approximations. 



5.4 Higher order 

The systematic expansion of the thermodynamic potential ip in powers of A can be carried 
out easily to higher orders. However the result involves some more detailed properties of the 
liquid at the effective temperature T*. For instance at second order one needs to know not 
only the free energy and pair-correlation of the liquid at temperature T*, but also the three 
points correlation. It is certainly interesting to try to push this expansion further, taking the 
information on the liquid at temperature T* from some numerical simulations. In this paper we 
have decided to stay within some relatively simple schemes which require only the knowledge 
of the pair-correlation g*(r). Therefore we shall not pursue this higher order expansion here, 
leaving it for future work. 



5.5 Harmonic resummation 

One can obtain a partial resummation of the small cage expansion described above by integrat- 
ing exactly over the relative vibration modes of the molecules. We shall use such a procedure 
here, which is a kind of harmonic expansion in the solid phase. 

We work directly with 1/a — and start from the replicated partition function (|3"5|), within 
the quadratic expansion of the interaction potential v in the relative coordinates u^. (Clearly it 
is assumed that the 1/a — >• + limit has been taken, and that its effect is to build up molecular 
bound states). The exact integration over the gaussian relative variables gives: 



m 



Nd/2 v /27T Ard(m ~ 1) f N 



Y[ d d Ti exp —/3m ^2 v ( r i ~ r j) 



m 



1 



Tr log {(3M) 



(41) 
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where the matrix M, of dimension Nd x Nd, is given by: 



M (in)Uf) = hi v ^i. r i ~ r k) ~ V^in - Tj) (42) 
k 

and v^ v (r) = d 2 v/dr fl dr l/ . We have thus found an effective Hamiltonian for the centers of 
masses r« of the molecules, which basically looks like the original problem at the effective 
temperature T* = T/m, complicated by the contribution of vibration modes which give the 
'Trace Log' term. We expect that this should be a rather good approximation for the glass 
phase. Unfortunately, even within this approximation, we have not been able to compute the 
partition function exactly. The density of eigenstates of the matrix M is a rather complicated 
object and we have developed a simple approximation scheme in order to estimate it. 

We thus proceed by using a 'quenched approximation', i.e. neglecting the feedback of 
vibration modes onto the centers of masses. This approximation becomes exact close to the 
Kauzman temperature where m — > 1. The free energy is then: 

d dim — 1 I I m — 1 

m™,T) = -— log(m) - ^log(27r) - l og Z(T*) + — — (Tr log (/?M)>* (43) 

2m 2m mJM 2m 

which involves again the free energy and correlations of the liquid at the temperature T*. 
Computing the spectrum of M is an interesting problem of random matrix theory, in a subtle 
case where the matrix elements are correlated. Some efforts have been devoted to such compu- 
tations in the liquid phase where the eigenmodes are called instantaneous normal modes f49j. 



It might be possible to extend these approaches to our case. Here we shall rather propose a 
simple resummation scheme which should be reasonable at high densities-low temperatures. 

Considering first the diagonal elements of M, we notice that in this high density regime 
there are many neighbours to each point, and thus a good approximation is to neglect the 
fluctuations of these diagonal terms and substitute them by their average value. We thus 
write: 

1 f 

J2 v ^( r i ~ r k) — <V j / d d rg*(r)Av(r) = r (44) 
k d J 

Here and in what follows, we have not written explicitly the density. We choose to work with 
density unity in order to simplify the formulae: this value can always be obtained by using an 
appropriate scale of length. In the approximation ([44]) the diagonal matrix elements are all 
equal and can be factorized, leading to: 

(Trlog(f3M))* = Ndlog((3r ) + (Tr log - ^(r* - r fc ))>* (45) 

This form lends itself to a perturbative expansion in powers of 1/tq. The computation of the 
p-th order term in this expansion, 

T p = (-l) p_1 overpr p ( J2 "wwfci ~ r i 2 )--V-iM P ( r Vi ~ ^K^iK, - r h) ( 46 ) 

h--*p 

Ml ■••Mp 

still involves the p-th order correlation functions of the liquid at T*. We have approximated 
this full correlation by introducing a simple 'chain' approximation involving only the pair 
correlation. This chain approximation consists in replacing, for p > 2, the full correlation by 
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a product of pair correlations. It selects those contributions which survive in the high density 
limit; systematic corrections could probably be computed in the framework of the approach of 
[]5U| ], we leave this for future work. Within the chain approximation, T p is approximated by: 



J dx 1 ...dx p g*(xi, ....,x p ) [v ^{xi - x 2 )...v flp _ llJ , p (x p - 1 - x p )v^ vill {x p - xi)\ 

'J,p 

J dx 1 ...dx p [g*{xi - x 2 )v fillX2 (x 1 - x 2 )]...[g*(x p - x 1 )v IJvlll (x p - xi)] . (47) 



iii... tip 



iii-.-np 



In this last form we need to compute a convolution which can be factorised through the 
introduction of the Fourier transform of the pair correlation function. We thus introduce the 
Fourier transformed functions a and b which are defined from the pair correlation g*(r) by: 

d d r g*(r)v, u (r)e ikr = 5, u a(k) + - -^j b(k) . (48) 

In terms of these Fourier transforms, the p-th order term in the l/r expansion is simply 

_ r d d k ( .,, rf-l 7/ ,A P s r d d k ( .,, 1 7/1 A p . , 

T ' = I W¥ [ a{k) + — m J + {d - l) I W? \ a(k) - 1 m ) • m 

and the summation of the series over p is easily done, so that the free energy per particle within 
the chain approximation of the harmonic resummation is: 

/30(m, T) = -± log(m) - log(2vr) - log Z(T*) + log(/?r ) 

+ (m-1) r d d k /J g(k) + ^b(k) \ | (a{k) - \b{k) \\ 



2m J (2tc) 3 \ \ r J \ r o J J 

(m-1) f d „^(r) 2 
- ^- / d d rg(r) > - \ (50) 

where the function L3 is defined as: 

L 3 (x) = log(l -x) +x + x 2 /2 (51) 

We can thus compute the replicated free energy F m only from the knowledge of the free energy 
and the pair correlation of the liquid at the effective temperature T*. The results will be 
discussed in section |7|. 



6 A systematic Approach: molecular HNC closure 
6.1 Density functional 

As we have seen before, one can choose as an order parameter the generalised inter-replica 
correlation, deduced from the original partition function by the functional derivative: 

oir 1 r m ) - 1 6hgZm (52) 
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In order to study the free energy at fixed order parameter, one can perform the functional 
Legendre transform: 

Mp] = -— logZ m [0] - — / dr l ...dr m p(r l ,...,r m )W(r\...,r m ) (53) 
m m J 

and the aim is to optimize this new function with respect to p. 

In the ideal case where there are no interactions, this thermodynamic potential is: 

T r ntr 1 r m ) 



^ d [p] = — / dr\..dr m p(r\ r m ) log rK - (54) 

m j e 

We need to add to this piece the part which comes from the interactions. This is non trivial; 
in the next section we shall use the HNC approximation for this function. 

6.2 Molecular HNC equations 

The free energy in the HNC approximation is derived in the appendix I. It is a functional of 
the molecular density p(x) and the two point correlation g(x,y) = 1 + h(x,y). Here and in 
the following, the letters x, y and z without any index denote md— dimensional vectors (e.g.: 
). The molecular density is our order parameter. The result for ijj is: 

Pfp = y~ J dxd yp( x )p(y) ld( x ' y) lo &9(x, y) - g(x, y) + l + 0v(x, y)g{x, y)} 

If 1 \ 1 r p(x) 

— Tr log(l + hp) — hp A — hphp I H / dxp(x) log (55) 

2m V 2 / m J e 

where the potential is v(x,y) = J2 a v i x a —y a )- In the trace term all products are convolutions. 
For instance the lowest order term in the small p expansion of the trace is: 



d md x d md y d md z h(x, y)p(y)h(y, z)p(z)h(z, x)p(x) (56) 



We would like to optimize the thermodynamic potential tp with respect to the molecular 
density p(x) and the two point function g(x, y). We shall work at low temperatures for which 
p should be nearly gaussian. We thus choose an Ansatz for p of the type (always with a choice 
of average density equal to one): 



p(x) = j*xn(^ { - { *z£' ,i2A)) ) 



a=l 

/2nA\ d / 2 



\ m 



) 2{2nArmd/2eXV {-^^ Xa - xb)2 ) 



(57) 



where the molecular density is parametrized by the single parameter A. 
The ideal gas contribution (last term in (p5|) gives: 

/ [] d d x a p(x) log — = N [ ^(1 — m) log(27rA) + ^(1 - m) - ^ logm - 1 j (58) 
a e \2 2 2 J 

The interaction term is more complicated, and we have only succeeded in optimising it in 
the small cage regime. 
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6.3 Second order small cage expansion 

Here we shall solve in general for g in the limit of small cage radius A, expanding in powers of 
A. 

As usual we go to molecular coordinates, introducing x a = X + u a and y a = Y + v a , with 
the constraints: J2 a ua = J2 a v<x = 0- The molecular density fl57D depends only on the relative 
coordinates: 



P {u) ^ Po m^(E« a ) (^) d/ \2nA)-^e W (-^EK- 6 ) 2 ) 



(59) 



The w's are thus gaussian distributed with a second moment: 



(«> = A(5 ab --)5, u (60) 



m 



We shall expand the two point correlation in powers of the relative coordinates, using the 
notations: 

g ({x + u a },{Y + v a }) = g(x-y) + y, s ^( x - y )(t,K u : + v X]- 2K ^) 

+ E ^{X - Y) (yXK - <)K - <)] - 21f J (61) 
where the constant is chosen in such a way that, for any A: 

J dup(u) J dvp{v) g(X + u\ X + u m ; Y + v\ Y + v m ) = G(X - Y) (62) 
The constant turns out to be: 

= A(m - 1) V . (63) 

It is not difficult to see that, thanks to the constraint (^), the knowledge of the functions 
S and T is enough to compute the free energy to order A 2 . This computation is done in the 
appendix II. Here we just give the result. We write the free energy to second order in the form: 

^ = (3F + (3F^ + + (3F 2 (64) 

The zeroth order terms are: 

dl — m , , >N id — 2)1 — m d , . . 

(3F = log (2ttA) + K —— -L — logm 65 

2 m 2 m 2m 

PK = i/(0(- 1 °g( 1 + ^)) + ^)-^ A; ) 2 / 2 

1 r*. 



+ / d d r (G(r) log G(r) - G(r) + 1 + f3mv(r)G(r)) (66) 

2m J 

where H(r) = G(r) — 1, and H(k) is the Fourier transform of H(r). It is clear from(|66|) that 
the zeroth order correlation function G(r) is exactly the pair correlation of the liquid at the 
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effective temperature T* = T/m in the HNC approximation, we thus recover our previous 
results. 

The first order correction is: 

&F X = PA^f^ I d d rG(r) £ v^{r) (67) 

At this order we can easily optimize the free energy with respect to G(r), and with respect to 
the cage size A. We get back the same result for A and the free energy as we had in the direct 
first order small cage expansion (fH)]). 

The advantage of this molecular HNC approach is that we can compute the second order 
term without needing to solve for three point correlations in the liquid. The second order 
correction is: 

f3F 2 = A 2 ^/^r-^5:(^(r) 2 + 2^(r)T^(r) + 2r Ml/ (r) 2 ) 
+ A 2 ^± [d d rJ2(S,u(r) + 2T flu (r))(3v, u (r) 



fJ,V 

I'll 



The stationarity conditions on S and T are easily solved. One finds: 

T, v = - l -G(r)f3v, v (r) (69) 
while S + T is the solution of the linear equation: 

+ >- = 1 S? e ' k ' {S - {k) + T - (k)) Trm (70) 

The equation for G is also easily found. Expanding G = Go + AG\, one sees that Go is the 
pair correlation g* of the liquid at temperature T/m, while the correction G\ is the solution of 
the linear equation: 



Gi(r) , . v r d*k ikr H (k)(2 + H (k )) 

Go{r) +m I) ^v^r) - J {1 + Ho{k) y 



Gx(k) (71) 



The solution of these equations and the physical consequences are discussed in the next 
section. 

7 Results 

In this section we indicate how to obtain the thermodynamic properties of the glass within 
each of the previous approximation scheme, and we give the results. 
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7.1 Methodology 

We have developed in this paper three approximation schemes. 

The small cage expansion has been carried out directly to first order in section |5.3|, and 
agree with the first order expansion within the molecular HNC approach. Within this first 
order approximation, the cage size is given explicitly in ( |40|) and the corresponding free energy 
<j)(m) is given in (|39|). We need to study the m dependance of (p. Clearly the only ingredients 
we need are the free energy and pair correlations of the liquid at the temperature T* = T/m, 
which is a temperature which lies in the range of the glass transition temperature, as we shall 
see. These properties of the liquid could be obtained by various means; here we have used the 
HNC closure for the pair correlation and the corresponding free energy in order to get them, 
(obviously one could try to use better schemes of approximation for the liquid, depending on 
the form of v (r), in order to improve the results; our point here is not to try to get the most 
precise results, but to show the feasibility of a quantitative computation of glass properties 
using the simplest approximations). Given the temperature T, the procedure is the following: 
we vary the value of m, and for each value we can compute the cage size A and the free energy 
0(m). As expected on general grounds (see section [|), we find a free energy which increases 
with m until it reaches the critical value m*(T) (such that ( |TTD holds), which is the phase 
transition boundary. This critical value is defined by dip /dm = 0. The configurational entropy 
is given by the solution of the two general equations (0), and the free energy of the glass is 
nothing but <p(m*,T). We get the internal energy and specific heat by differentiating the free 
energy. The critical (Kauzman) temperature Tk is defined by m*(Tx) = 1. 

The second approximation scheme is the harmonic resummation method. Again we have 
an explicit form ([5(]) for the free energy per particle 0(m) only from the knowledge of the free 
energy and the pair correlation of the liquid at T* . Having this m dependance the procedure 
to get the thermodynamic results is entirely the same as that of the first order result. 

The third approximation scheme is obtained by the expansion of the molecular HNC free 
energy to second order in the cage size, as described in section ||. For given values of the 
temperature T and the number of replicas m, we first solve the standard HNC equations 
giving the pair correlation G(r) = g*(r) at the temperature T* = T/m. Then we can compute 
the functions S, T and the correction to the correlation G\ by solving the set of linear equations 



( |69| , |70| , |7TD . The free energy is then computed to second order as in (|64|). 

We use the results of the second order term in the expansion in a perturbative way which 
we shall now describe. One might be tempted to use the free energy computed to order A 2 
without expanding the solution to order A 2 . However this procedure is not useful because the 
equations truncated at the order A 2 do not have a solution. One must do the computation 
fully perturbatively in a consistent way, which we now explain. Let us define the various terms 
in this free energy as 

M(A, m) = 7o + Ay x + A 2 l2 + 73 log A (72) 
where the 7's are functions of m that we can compute. We suppose that the 72 term is small 
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and write the value A max which maximises [] the free energy as: 

2 

A - — — — ol3ll (7^\ 

■"-max % \ 1 U J 

7i 7i 

giving a free energy on this maximum approximately equal to 

ip(A max , m) = ^i(m) + V2M (74) 

with 

^i(m) = 7o-73 + 73log(-73/7l) 

^a(m) = 7273/71 (75) 

where ip2 is the correction term. This is a function of m which we maximise in order to find 
the critical value m*. Writing m* = m\ + 1712, where mi is the critical value computed to first 
order and m 2 is the correction, these numbers satisfy the equations: 

am 

m 2 = (^(m!)] (76) 

am y <7m z y 

For consistency of this perturbative expansion, one should then compute the saddle point value 
of A as: 

A _ 73(^1) 2 72(^i)73(^i) 2 m ^ d 73 (mi) 
71 (mi) " 7i( m i) 3 (9mi7i(mi) 
and the free energy of the glass as: 

if) = ^i(mi) +ip2{mi) (78) 

Having the free energy as a function of m we proceed as before by maximising it, following 
exactly the same steps as for the first order computation. 

7.2 Numerical procedure 

We have studied the case of soft spheres in three dimensions interacting through a potential 
v{r) = 1/r 12 . We work for instance at unit density, since the only relevant parameter is the 
usual combination T = pT^ 1 ^ 4 . 

For each of the three approximation schemes mentioned above, we need to compute the 
free energy and the pair correlation of the liquid in a temperature range close to the glass 
transition. We have used the HNC approximation to get both g(r) and the free energy. We 

^^One must maximise the free energy with respect to A, instead of the usual minimization procedure, when- 
ever m is less than 1. This is a usual aspect of the replica method, which is here a consequence of the fact that 
the free energy is proportional to to — 1 
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have solved the HNC closure equations numerically. For spherically symmetric functions in 
dimension three we use the Fourier transform for the radial dependance, in the following form: 



We discretize this formula introducing in r space a cutoff R and a mesh size a. In this 
way we have a simple formula for the inverse Fourier transform and we can also use the fast 
Fourier transform algorithm. In most of the computations we have taken a = 1/32.5 and 
L = 128 *asi4. We have checked that dividing a by 2 and multiplying L by two (thus going 
up to 512 points) does not alter the results. The solution of the equations can be found either 
by using a library minimization program , or a program which solves non linear equations. We 
have found first the solution at low enough density and then followed it by continuity while 
gradually increasing the density. 

The second order expansion of the molecular HNC theory requires some more work, because 
we need to compute the various tensors S^, T^, and the correction to G. After decomposing 
the tensors in their various irreducible components, using rotation invariance, these components 
are discretized on the same grid as g(r) and the linear equations are solved by a standard library 
routine. 

7.3 Critical temperature and effective temperature 

We plot in fig. f| the inverse of the effective temperature T*, equal to m*/T, versus the 
temperature T of the thermostat. The transition temperature is given by T* = T. This gives 
the ideal glass transition temperature. Within the first order expansion we find Tk — -14; the 
harmonic resummation gives Tk ~ .19 and the second order perturbation theory is Tk ~ .18 
We see that the two best methods, the second order and harmonic resummation, are in good 
agreement and give a critical value of Y around V ~ 1.52. This value of V is in good agreement 
with the published values of the glass transition of the soft sphere system, which range around 



We also notice that the effective temperature stays relatively constant when the actual 
temperature varies. Our results are not so far from a situation in which one would have 
T* ~ Tk, independently from the value of the temperature T, which means that m ~ T /Tk- 
A nearly linear variation of m versus T is often found in discontinuous spin glasses, where it is 
characteristic of a free energy landscape which is totally frozen in the whole low temperature 
phase 0]. It is worth noticing that such a relation has also been found for the temperature 
dependance of the fluctuation dissipation ratio (although, as this ratio is a dynamical quantity, 
it rather equals T/Tp, where To is the dynamical (mode-coupling) transition temperature). 

7.4 Cage size 

In replica space the cage size characterizes the size of the molecular bound state, in the approx- 
imation of quadratic fluctuations, as defined in (|26|). Its physical meaning is easily established: 
In the glass phase at low temperatures one can approximate the movement of each atom as 
some vibrations in a harmonic potential in the neighborhood of a local minimum of the energy. 
The typical square size of the displacement is given by: 




(79) 



1.6 0. 



A = ((r, - (r,)) 2 ) 




26 



0.20 



0.15 ■ 



§ 0.10 



0.05 ■ 



0.00 ' ■ " ■ ' ' " ' 1 

0.00 0.05 0.10 0.15 0.20 

T 

Figure 4: The effective temperature of the molecular liquid at the transition, T* = T/m*, 
versus the temperature T, computed in an expansion to first order (dashed-dotted line) and 
second order (full line) in the cage size A, and in the harmonic resummation (dashed line). 



which is the physical definition of the square size. The cage size is plotted versus temperature in 
fig. [5]. The cage size is nearly linear in temperature, as it would be in a T-independent quadratic 
confining potential. This indicates that the local confining potential has little dependance on 
the temperature in the whole low temperature phase. 



7.5 Free energy and specific heat 

In fig. P we plot the free energy versus the temperature for each of our three approximations. 
The strong consistency of the second order small cage expansion and the harmonic resummation 
are clearly seen. The data extrapolates at zero temperature to a ground state energy of order 
1.95. This is related to the typical energy of the amorphous packings of soft spheres. More 
precisely, if we consider all the amorphous packings of soft spheres at unit density, we can 
count them through the zero temperature configurational entropy. The lowest energy at which 
one can find an exponentially large number of such packings is the ground state energy of the 
glass state which we find within our approximations equal to 1.95. This could be amenable 



to some numerical test J5^, |53], |54]j . However in order to do such a test one must remember 
that we have not taken into account the existence of a crystal: therefore one must first remove 
all crystal like solutions, i.e. solutions which correspond to a crystal with some local defects. 
These solutions can be characterized by the presence of delta functions at the appropriate 
values of the momenta. This procedure of identifying crystal like solutions has been explicitly 



done numerically in | 54fl . Generalizing the present result to hard spheres would allow for a 



computation of random close packing density, a notion which is often used in granular materials 
In fig.^ we plot the internal energy of the glass versus temperature, computed in each of our 
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Figure 5: The parameter A/T versus the temperature, computed in an expansion to first 
order (dashed-dotted line) and second order(full line) in the cage size A, and in the harmonic 
resummation (dashed line) 




Figure 6: The free energy versus the temperature, computed in an expansion to first order 
(dashed-dotted line) and second order(full line) in the cage size A, and in the harmonic resum- 
mation (dashed line). 
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Figure 7: The internal energy versus the temperature, computed in an expansion to first 
order (dashed-dotted line) and second order(full line) in the cage size A, and in the harmonic 
resummation (dashed line). Also shown is the internal energy of the liquid (dotted line). 

approximation schemes. Also shown is the internal energy of the liquid. The internal energy 
is continuous at the transition. 

In fig. H we plot the specific heat versus temperature. It is basically constant and equal to 
3/2. The fluctuations are numerical errors due to the extraction of the specific heat through 
the numerical second derivative of the free energy. A specific heat C = 3/2 is nothing but the 
Dulong- Petit law (we have not included the kinetic energy of the particles, which would give an 
extra contribution of 3/2). This result is very welcome: in fact if we had treated the crystal at 
the same level of approximation as we considered here for the glass, we would get the Einstein 
model for which the specific heat is also given by the Dulong-Petit law. Thus we have found 
that the specific heat of the glass is equal to that of the crystal, which is a good approximation 
of the existing data. Notice that it was not obvious at all a priori that we would be able to get 
such a result form our computations, since we are performing some computations purely in the 
liquid phase, with a liquid pair-correlation etc... The fact of finding the Dulong-Petit law is an 
indication that our whole scheme of computation gives reasonable results for a solid phase. At 
a later stage we would like to go beyond the Dulong-Petit law and get a better computation 
of the spectrum of soft vibration modes in order to get a Debye-like law. This is left for future 
work. 

7.6 Configurational entropy 

In fig. |S] we show the configurational entropy versus the free energy at various temperatures, 
including the zero temperature case. We have included here for simplicity only the result from 
the harmonic resummation procedure. 

We notice that the various curves corresponding to different temperatures are not far from 
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Figure 8: The specific heat of the glass versus the temperature, computed in an expansion 
to first order (dashed-dotted line) and second order(full line) in the cage size A, and in the 
harmonic resummation (dashed line). The dotted line is the specific heat of the liquid. 
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Figure 9: The configurational entropy £(/) versus the free energy, computed within the har- 
monic resummation, at temperatures T = 0., .05, .1 (from left to right). 



30 



being just shifted one from another by adding a constant to the free energy. This indicates 
that the main effect of temperature amounts to an additive constant in the energies of all 
amorphous packings. This would be the case if the states at finite temperature could be deduced 
continuously from the zero temperature amorphous packings, with an extra contribution to 
the free energy coming from the vibrations, if the vibration spectrum is more or less state 
independent. 

7.7 Dynamical transition 

As we discussed in the introduction, at the mean field level there exists a dynamical transition 
at a temperature T D larger than the thermodynamic transition temperature T K . This phase is 
characterized by the dynamic statement that a system will remain forever in the same valley, 
and its free energy is greater than the equilibrium one because it misses the contribution of the 
configurational entropy. It is thus evident that this dynamic phase is just a mean field concept, 
which should disappear when corrections, such as activated processes, due to the short range 
nature of the potential, are taken into account. However if the barriers are sufficiently high, 
metastable states have a very large life time and they strongly affect the dynamics. It would 
be thus interesting to try to compute the 'dynamic transition temperature' in these systems. 

In the framework of the harmonic resummation one finds that the approximation breaks 
down at small but positive e if the matrix of second derivatives has negative eigenvalues. 
From this point of view the appearance of negative eigenvalues signal the dynamic transition. 
Unfortunately in our chain approximation all the eigenvalues are positive at all temperatures 
and no dynamic phase transition can be seen: the free energy is always well defined for small e. 
This negative result is due to the fact that the chain approximation we use may be reasonable 
at low temperature but it is certainly not good at high temperatures. This problem will 
disappear if one uses a better method to compute the spectrum, giving reasonable results also 
at higher temperatures. On the other hand in the framework of the small cage expansion the 
perturbative method assumes that there is always a bound state. Although this should not be 
true at high temperature, the breakdown of this assumption cannot be seen in a perturbative 
approach. 

It is clear that a study of the dynamical phase transition should be done using some 
different tools than the one we have developed here. This is not surprising: the dynamical 
phase transition is present at a temperature higher than the static one and the approximations 
which we have been using are low temperature ones. 

8 Discussion and perspectives 

Deducing the thermodynamic properties of the glass from those of a liquid may look crazy. Of 
course the main trick is that we use a molecular liquid, with a variable number m of atoms per 
molecule, which will have a glass transition at a temperature lower than Tk whenever m < 1. 
We wish to underline again what is the basic hypothesis of our approach. We assume that there 
exists a thermodynamic glass transition, which is of the general type described in our 'basic 
scenario'. This assumption means that there exists a path in the m , T space which connects 
the points m < 1 , to the high temperature region without crossing any transition. 
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If this is true (and this is known to happen in many models) the situation is rather simple 
and corresponds to what is called in the literature one step replica symmetry breaking. This 
situation corresponds to the case in which the deep minima of the free energy are completely 
uncorrelated |22], f46|. One could think of checking this hypothesis numerically by computing 
for small systems all the metastable states at zero temperature, and studying the distribution 
of their energies. Let us mention for completeness that there exist models in which the deep 
minima of the free energy are partially correlated (this is very probably the case of spin glasses 
47], |4~8|). In such a case any path in the m — T space which connects the point (m, T"( m )) to the 



high temperature region crosses a phase transition, and one would need to introduce a more 
complex construction in order to avoid this singularity 

The approach described in this paper opens the way to the computation of the thermody- 
namic properties of glasses at all temperatures using the generalization of the standard tools 
of liquid theory. Although it is not explicitly discussed in this paper, this approach allows 
also the computation of the density correlation function g(r) in the glassy phase; we plan to 
address this point in the next future. 

It is clear that the results presented here just use the simplest possible non trivial approx- 
imations. Nevertheless, within these simple approximations, we have shown that a reasonable 
value of the Kauzman temperature can be derived, as well as several thermodynamic properties 
of the glass phase: the internal energy, free energy, configurational entropy and specific heat, 
and the cage radius. Obviously our study so far has been restricted to equilibrium proper- 
ties, and the equilibrium situation is very difficult to reach experimentally. However one can 
think of measuring each of the above properties in numerical simulations, where the joint use 
of smart algorithms and small enough system can allow to thermalize. The extension of the 
present methods to binary mixtures is a work that must be done in order to allow for a more 
precise comparison with the results of numerical simulations. Some steps have already been 
done in this direction |40] . 



This equilibrium study is to be considered as a first step before dealing with the out of 
equilibrium dynamics. Beside the dynamics in the low temperature phase, a very interesting 
and open problem is the computation of the time dependent correlation functions (and as a 
by-product the viscosity) in the region above Tk- However a better understanding of activated 
processes in this framework is a crucial prerequisite. 

Within the equilibrium framework, we have implemented so far our general strategy using 
rather crude methods. These should be improved, which means that one must perform a more 
careful study of the molecular liquid. There are many directions in which one could move: 

• Improve the computation of the spectrum in the harmonic approximation. This harmonic 
approximation should be excellent and allow to study from first principles all the low 
temperature anomalies which have been observed in glasses. Within this approximation 
one just needs to study the liquid of the centers of masses of the molecules, which interact 
through the effective interaction described in (pLT|). Of course the interaction term coming 
from the Tr log term is not easy to deal with, but still this is a very well defined problem 
of liquid theory for which precise approximation scheme should be developed. 

• Use approximations different from HNC, which may work better in the liquid phase. 
Obviously this will depend on the interaction potential, and a detailed study of several 
different types of potentials would be very interesting. 
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• Use numerical simulation in the liquid phase in order to get some higher order coefficients 
of the A expansion: these are given by higher order correlation functions which could be 
measured in simulations. 

• Introduce resummation techniques that are more efficient than the harmonic one. 

Some of the previous described techniques could also be used to understand better the prop- 
erties of the dynamical phase transition. 

To summarize, our approach transforms the problem of the thermodynamics of the glass 
phase into a problem of a (complicated) liquid state. We hope that the sophisticated methods 
developed in liquid state theory will be brought to bear on the study of glasses. 
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Appendix I: HNC closure 

For completeness, we give here a derivation of the HNC free energy (|53|) for our molecular 
replicated system. One could use the standard diagrammatic method p6[ , but here we shall 
follow the 'cavity' like method of Percus [57|. We study N molecules with coordinates Xi,i G 
{l,..,iV}. Each Xi stands for the coordinates of all atoms in molecule i: Xi = {x"},a G 
{1, m}. The energy of the system is given by 

E = Y.V{x i ,x j ) + Y. u ( x i) ( 81 ) 

i<j i 

where v is the intermolecular potential (in our case we would have V(x, y) = J2a v i xa ~ U a ) 
but we shall keep a general V in this appendix), and the external potential u(r) has been 
introduced for future use. 

We shall need the following definitions. The one molecule density is 

p(*)=E<n^-* a )> , (82) 

i a 

where the average (.) is with respect to the Boltzmann measure exp(— f3E). The two molecules 
correlation is: 

P v (x, y) = £<n K< - * a ) II -* b )) = p(x)g(x, y)p(v) m 

i+3 a b 

where we have also defined the pair correlation function g(x,y), which goes to one at large 
(center of mass) distance. The connected pair correlation is: 

h(x,y)=g(x,y)-l (84) 
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Elementary functional differentiation gives 

dp(x) 



p(x)5(x-y) + p(x)h(x,y)p(y) (85) 



d(-Mv)) 

One can also introduce the direct correlation function c(x, y) through: 

d{-j3u{x)) 1 



dp(y) p(x) 



S(x -y) - c(x,y) (86) 



the direct correlation is thus related to the connected pair correlation through the Ornstein- 
Zernike equation c = (1 + hp)~ 1 h which reads more explicitly: 



c(x,y) = h(x,y) + J dx\h(x, x\)p(xi)h(xi, y) + / dxidxzhfa, x\)p(xi)h(xi, X2)h(x2, y) + ... 

(87) 

The idea of Percus is to compute the pair correlation by considering the one point density 
with a molecule fixed at one point. Let us consider a problem in which we have added one extra 
molecule, fixed at a point z = {z 1 , ...,z m }. This extra molecule creates an external potential 
u(x) = V(x,z). The one point density in the presence of this external potential, p u (x), is 
related to the density p(x) and pair correlation g(x,z) in the absence of an external potential 
through the conditional probability equation: 

,( 2 ) 



Pu(x) = p ( >(x, z)/p(z) = p(x)g(x, z) 

In order to try to build a successful approximation scheme, let us introduce two quantities 
R u {x) and S u (x) which we can calculate in presence of the external potential, or when this 
potential is turned off (u = 0). If their variations are smooth enough, one can approximate 
their variations by the first order term: 

R u (x) ~ R u=0 (x) + I dy 5 -^\ (S u (y) - S u=0 (y)) (89) 

The standard perturbation theory would be obtained by taking R u (x) = p u (x) and S u (x) = 
u(x). However the linear truncation (1891) can be better behaved with some better choices of 



the functions R and S. The HNC closure corresponds to taking f57 |: 

R u (x) = log (p u (x)e^) ; S u (x) = p u (x) . (90) 

Then we have 

^|(n = 0)=c(x,y) (91) 
and the linear equation (|89|) becomes: 

\ogg(x,z)+pV{x,z)= dyc(x,y)p(y)h(y,z) (92) 



Together with the inversion relation (p7|), this defines a closed set of equations for the one 
and two point molecular densities which are the HNC closure. It is easy to show that these 
equations express the stationarity of the free energy functional i[>[p,g] defined in d55|) , with 
respect to variations of g. 

The result for the free energy can be deduced if we assume that: 
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• There exists a variational principle where the free energy is a functional of g and p. 

• The potential (3V(x) enters in the free energy is such a way that the internal energy takes 
the exact form 1/2 / dydxp(x)p(y)g(x,y)V(x,y). 

• The free energy functional at g — 1 and v — 0, which depends only on p is given by the 
exact form 

— / TT d d x a p(x) log (93) 
m J a e 

These three conditions fix in a unique way the free energy functional and are satisfied in the 
previous approach. Indeed the second condition implies that the free energy ip can be written 

as 

pip = (3/2 / dydxp(x)p(y)g(x, y)V(x, y) + x[g, p] (94) 



where x does not depends explicitly on (3. If we differentiate the previous equation with respect 
to g we find 

(3/2p(x)p(y)V(x,y) + -^—. (95) 

og{x,y) 

If we identify the previous equation with eq. (|55|) (after multiplication by p(x)p(y)) we find that 
the proposed free energy (eq. (|9"2"D) has the same derivative with respect to g of the exact one. 
Now the only ambiguity that remains in the free energy is its value at g — 1 and v = 0, which 
is fixed from the condition (3). 



Appendix II: Second order small cage expansion 

Here we carry out the small cage expansion of the molecular HNC equations to second order. 
We start from the HNC free energy (|55|), we introduce the center of mass and relative coordi- 
nates, x a = X + u a and y a = Y + u a , and we expand in the cage size A, using the molecular 
density (59) and the decomposition of the correlation function given in (|6lD . 



We shall examine successively the various pieces of 2m(3ip. The form of the simplest piece 
is deduced trivially from the constraints (|62j) : 



J dxdy p(x)p(y)(l-g(x,y)) = J dXdY (1 - G(X, Y)) (96) 

(we remind that here x and y stand for all the molecular coordinates and are therefore md- 
dimensional vectors, while the center of mass coordinates X and Y are <i-dimensional) . 
We go next to the piece involving the potential: 

W dxdy p(x)p(y) v(x a - y a )g(x, y) . (97) 

a J 

We expand the potential as: 

£ v(x a - y a ) = mv(X -Y) + W v, u (X - Y) - «•)(«; - «•) + ... (98) 
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and expand the correlation according to (|H]). Thanks to the constraint (0), the term in 
mv(X — Y) contributes exactly as: 

m J dXdY G(X - Y)v(X - Y) (99) 

to all orders in u, v. The term in v IMl/ (X — Y) contributes a piece of order A which is: 

f dXdY p(u)dup(v)dv G(X - Y)J2v„AX - Y) " <)« - <) (100) 

p,v a 

and a piece of order A 2 which is: 

f dXdYp(u)dup(v)dvJ2v^(X ~ Y) EK - - <) ( 101 ) 

The last piece of order A 2 comes from the fourth derivative of v in (|98|) : 

dxdYp(u)du P (v)dvG(x-Y) J2 v^(^-^EK-<)«-<)K-<)K-<) ( 102 ) 



/ivpcr 



Notice that the use of (|62D allow us to find the order A 2 expression without ever introducing 
the order A 2 term in the expansion of the pair correlation. This will also be true for the 
other contributions below. This strategy is crucial for keeping the computation not too big. 
The various pieces are now easily computed using the fact that u and v variables are gaussian 
distributed with the second moment given in (|60"1) . We get: 



PY,J dxdyp(x)p(y) v(x a - y a )g(x,y) = V J dXG(X) (mv(X) + A(m - l)J2 v w( X )) 

+VA 2 (m-l) J dXj(2S p (%(X) + ff F (%m + ^ ww (Z))l03) 

We now turn to the 'glogg' term in the free energy 2m(3ip. Expanding as before, we get: 
J dxdyp(x)p(y)g(x,y) log g(x,y) = J dXdYG(X - Y) \ogG{X - Y) + (104) 
+ f dXdY p(u)dup{v)dv — — \ — — £ 

S,u{X - Y) + vy„] - 2K^j + T, U (X - Y) (fe« - W J)(«* - v b v )\ - 2K^ 

s p m - y) (x>x + «\ - 2R P °) + T ^ x - y ) (xk - v %< - o\ ~ 2 ^) 
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which gives after performing the gaussian u and v integrals: 
V J dX \G(X) logG(X) + 4A2(m ~ 1 



G{X) 



E 



-S^WS^X) + S^ V {X)T^{X) + T pp {X)T pu {X) 



We now study the last piece of 2m/3^>, namely the convolution term 
00 (— 1)p /■ 

/ dx\...dx p p(x\)h(xi, x 2 )p(x 2 )h(x 2 , x 3 )...p(x p )h(x p , x\) . 

p=3 P 



(105) 



(106) 



Here again each Xj is a m<i dimensional vector including all molecular coordinate, which we 
decompose into the center of mass Xj and the relative coordinates u®. Therefore each piece 
h(xj,Xj + i) in the above product is expanded as: 



h(xj, Xj + i) = h(Xj — X 



j'+i/ 



+ E 



V 6 



+T lxu (X j - X j+1 ) (]T[(< 
V b 



+i,„)] ~ 



(107) 



We notice again that higher order terms do not contribute to order A 2 . The second order 
terms generated by the expansion (|107|) when it is inserted into ( |106| ) are obtained by picking 
up the 'h(Xj — Xj+i)' contribution in all but two values of j. In order for the result not to 
vanish (because of (|62|) ) , we need that this two special values of j be neighbours. We thus get 
the following order A 2 contribution to the convolution term: 



00 . 

(— l) p / dX\...dX v p{u\)du\...p{u v )du v 

p=3 J 



E 



S^X, - X 2 ) (E[<„<* + u b ^u b 2 J - 2K,)j + 

T pv (X x - X 2 ) - - u\ v )] - 2K^ 

S pa {Xi - X 2 ) (EKXv + 4 tP u b J - 2K p ^ + 

T P Mi ~ X 2 ) (Y\{u\ jp - < p )« CT - u b >a )} - 2K p )j 
h(X 3 - X 4 )...h(X p ^ - X p )h(X p - X,) 



(108) 



After performing the gaussian u and v integrals, we find an expression in terms of the Fourier 
transformed functions h(k), S^k) and S pv {k): 



00 Jdh, 

2VA 2 {m-l)Y J {-l) P T, / j^Hky-'iS^+T^k))' 

P =3 \iv J y Z7T ) 



(109) 
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involving a simple geometric series. 

Grouping together all the pieces of the free energy tp which we have considered, we obtain 
the second order expression of the free energy used in ( |B4| - |5"8| ) . 
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